Epel recruits BET family bromodomain 
protein Bdf2 to establish heterochromatin 
boundaries 

Jiyong Wang/ Xavier Tadeo/ Haitong Hou/ Patricia G. Tu/ James Thompson,^ John R. Yates III/ 
and Songtao Jia^'^ 

^Department of Biological Sciences, Columbia University, New York, New York 10027, USA; ^Department of Chemical 
Physiology, The Scripps Research Institute, La Jolla, California 92037, USA 



Heterochromatin spreading leads to the silencing of genes within its path, and boundary elements have evolved to 
constrain such spreading. In fission yeast, heterochromatin at centromeres I and III is flanked by inverted repeats 
termed IRCs, which are required for proper boundary functions. However, the mechanisms by which IRCs prevent 
heterochromatin spreading are unknown. Here, we identified Bdf2, which is homologous to the mammalian 
bromodomain and extraterminal (BET) family double bromodomain proteins involved in diverse types of cancers, 
as a factor required for proper boundary function at IRCs, Bdf2 is enriched at IRCs through its interaction with the 
boundary protein Epel. The bromodomains of Bdf2 recognize acetylated histone H4 tails and antagonize Sir2- 
mediated deacetylation of histone H4K16. Furthermore, abolishing H4K16 acetylation {H4K16ac) with an H4K16R 
mutation promotes heterochromatin spreading, and mimicking H4K16ac by an H4K16Q mutation blocks 
heterochromatin spreading at IRCs, Our results thus illustrate a mechanism of establishing chromosome boundaries 
at specific sites through the recruitment of a factor that protects euchromatic histone modifications. They also 
reveal a previously unappreciated function of H4K16ac in cooperation with H3K9 methylation to regulate 
heterochromatin spreading. 
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In eukaryotes, genomic DNA is folded with histone and 
nonhistone proteins in the form of chromatin. Chromatin 
regulates diverse DNA-based processes such as gene 
transcription, recombination, and DNA damage repair 
(Bannister and Kouzarides 2011; Zentner and Henikoff 
2013). Based on the level of compaction, chromatin is 
categorized as euchromatin or heterochromatin. Euchro- 
matin is gene-rich and usually associated with active 
transcription, whereas heterochromatin is gene-poor and 
highly compacted and limits the access of transcription 
and recombination machinery (Grewal and Jia 2007; 
Beisel and Paro 2011). Heterochromatin has the intrinsic 
ability to spread in a sequence-independent manner, inac- 
tivating genes over long distances (Talbert and Henikoff 
2006; Grewal and Jia 2007; Moazed 2011). Thus, to 
maintain stable gene expression patterns, it is essential to 
protect euchromatin from encroachment by heterochro- 
matin (Gaszner and Felsenfeld 2006; Valenzuela and 
Kamakaka 2006). 
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Heterochromatin spreading depends on the self-propa- 
gation of heterochromatin-associated histone modifica- 
tions and chromatin proteins (Rusche et al. 2003; Talbert 
and Henikoff 2006; Moazed 2011). Histones within het- 
erochromatin are methylated on H3 Lys9 (H3K9me), 
which serves as a signal for the recruitment of hetero- 
chromatin protein 1 (HPl) family proteins (Rea et al. 
2000; Bannister et al. 2001; Nakayama et al. 2001). HPl 
recruits the H3K9 methyl transferase SUV39 (Schotta 
et al. 2002; Stewart et al. 2005), leading to methylation 
of adjacent nucleosomes. This positive feedback loop 
induces heterochromatin spreading outwards from its 
nucleation sites (Talbert and Henikoff 2006). Histones 
within heterochromatin regions are also hypoacetylated. 
In budding yeast, which lacks H3K9me and HPl proteins, 
the Sir2/Sir3/Sir4 complex mediates heterochromatin 
spreading in a similar fashion, where Sir2/Sir4-mediated 
deacetylation of histone H4K16 and recruitment of Sir3 
form a similar self-propagating mechanism for hetero- 
chromatin spreading (Kurdistani and Grunstein 2003; 
Rusche et al. 2003; Moazed 2011). Given that H4K16 
acetylation (H4K16ac) status directly affects the compac- 
tion levels of nucleosome arrays in vitro (Shogren-Knaak 
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et al. 2006), modulation of H4K16ac levels could poten- 
tially be a universal mechanism to regulate heterochro- 
matin assembly and spreading. However, functional ho- 
mologs of Sir3 and Sir4 are absent from organisms using 
H3K9me/HPl systems to assemble heterochromatin, and 
the role of H4K16ac in regulating heterochromatin spread- 
ing outside of budding yeast is unknown. 

Spreading of heterochromatin is limited by the avail- 
ability of silencing proteins and competition between 
positive and negative regulatory components. As a result, 
the heterochromatin-euchromatin border is often not 
precisely defined but instead exhibits an extended tran- 
sition zone (Kimura and Horikoshi 2004). In certain cases, 
heterochromatin is confined by specific DNA sequences 
termed boundary elements, which are characterized by 
sharp transitions in histone modification profiles (Gaszner 
and Felsenfeld 2006; Valenzuela and Kamakaka 2006). 
When the functions of such boundaries are disrupted, 
heterochromatin spreads and silences genes located out- 
side of boundaries. Thus, boundary elements are essential 
for maintaining stable gene expression patterns. The mech- 
anisms that specify heterochromatin boundaries are diverse, 
but most function by either recruiting histone-modifying 
enzymes to directly antagonize heterochromatin-associated 
histone modification cycles or tethering chromatin to the 
nuclear envelope to physically separate different chroma- 
tin domains (Gaszner and Felsenfeld 2006; Valenzuela and 
Kamakaka 2006). 

Heterochromatin assembly has been extensively studied 
in fission yeast. Constitutive heterochromatin is mainly 
present at the pericentric, subtelomeric, and silent mat- 
ing type regions of the fission yeast genome. All of these 
regions contain similar DNA elements that serve as 
nucleation centers of heterochromatin assembly through 
the RNAi pathway (Moazed 2009; Lejeune and Allshire 
2011; Castel and Martienssen 2013). Additional DNA- 
binding proteins such as Atfl/Pcrl and Tazl promote 
heterochromatin nucleation within the silent mating type 
region and telomeres, respectively (Jia et al. 2004; Kanoh 
et al. 2005). These nucleation sites recruit the histone 
methyltransferase Clr4 complex (CLRC), leading to H3K9me 
and the recruitment of the HPl homolog Swi6 (Nakayama 
et al. 2001; Hong et al. 2005; Horn et al. 2005; Jia et al. 2005). 
Once initiated, heterochromatin spreads from these nu- 
cleation centers into neighboring regions (Hall et al. 2002; 
Kanoh et al. 2005). The silent mating type and pericentric 
regions are surrounded by special DNA sequences that 
correspond to sharp transitions of histone modification 
profiles (Noma et al. 2001; Cam et al. 2005). At the mating 
type region, heterochromatin is flanked by two inverted 
repeats that recruit TFinC (Noma et al. 2006), the ma- 
chinery required for RNA polymerase III (Pol in)-mediated 
transcription (Paule and White 2000). TFIIIC tethers its 
binding sites to the nuclear periphery, a process that may 
contribute to formation of heterochromatin boundaries 
(Noma et al. 2006). Some pericentric heterochromatin 
borders correlate with the presence of tRNA gene clusters 
(Cam et al. 2005), which are critical for boundary function 
(Scott et al. 2006). Given that Pol HI transcribes tRNA 
genes, it is possible that these tRNA genes mediate 



boundary formation through a TFIIIC -based mechanism. 
There are also inverted repeats flanking centromeres I 
and III termed IRCs, which are also required for boundary 
activity (Noma et al. 2006), but the mechanism by which 
IRCs define heterochromatin boundaries is unknown. 

Epel is a JmjC domain-containing protein required for 
boundary function (Ayoub et al. 2003). It is recruited to all 
heterochromatin domains in a Swi6-dependent manner 
(Zofall and Grewal 2006; Trewick et al. 2007; Sadaie et al. 
2008) but is concentrated at borders of heterochromatin, 
such as IRCs, due to the degradation of Epel in the middle 
of heterochromatin by the proteasome, mediated by the 
Cul4-Ddbl E3 ubiquitin ligase complex (Braun et al. 201 1 ). 
Although JmjC domain proteins are generally involved in 
histone demethylation, Epel does not possess such activ- 
ity (Tsukada et al. 2006), and the mechanism by which 
Epel functions in blocking heterochromatin spreading is 
still unknown. 

By systematically screening a fission yeast deletion 
library for mutations that result in heterochromatin bound- 
ary defects, we identified a novel protein, Bdf2, required for 
heterochromatin boundary function at IRCs. Bdf2 interacts 
with Epel and is recruited to IRCs in an Epel -dependent 
manner. Through its double bromodomains, Bdf2 preferen- 
tially interacts with acetylated histone H4 tails and pro- 
tects them from Sir2-mediated deacetylation of H4K16. 
Our results demonstrate that protection of a euchromatic 
histone modification by a chromatin modification binding 
protein at a specific DNA element is sufficient to establish 
a heterochromatin boundary. We also uncovered a previ- 
ously unappreciated role of H4K16ac in regulating hetero- 
chromatin spreading in a system where heterochromatic 
silencing is dominated by H3K9me and HPl. 

Results 

IRCl is a heteTOchwmatin boundary 

The right border of pericentric heterochromatin on chro- 
mosome I corresponds to the presence of an IRC element 
[IRCIR). There are neither any tRNA genes nor any 
detectable TFIIIC enrichment (Noma et al. 2006), making 
it an ideal location to study IRC function. We inserted 
a ura4^ reporter gene to the right of IRCIR [IRCIR:: 
ura4'^) (Fig. lA). Because heterochromatin spreading de- 
pends on the dosage of heterochromatin proteins, over- 
expression of Swi6 has been used to improve heterochro- 
matin spreading at the silent mating type region and 
centromeres (Noma et al. 2001, 2006). Therefore, we 
inserted an additional copy of swi6'^ at the ade6 locus to 
enhance heterochromatin spreading. In wild-type cells, 
IRCIR ::ura4^ was fully expressed, and cells could not 
grow on counterselective medium containing 5-fluoroorotic 
acid (FOA), which is toxic to cells expressing Ura4 (Fig. 
IB). However, when IRCIR was deleted, the spreading of 
heterochromatin led to silencing of this reporter, as evi- 
denced by increased growth on FOA medium (Fig. IB). This 
silencing was eliminated in ch4A cells (Fig. IB), suggesting 
that Clr4-mediated H3K9me is critical for heterochroma- 
tin spreading. 
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Figure 1. Bdf2 is required for boundary function. {A,C,G) Schematic diagrams of the IRClR::ura4'^, L5-IRClR-ura4^, and IRC3L:: 
ura4'^ reporters. {B,D,H) Tenfold serial dilution analyses of the indicated yeast strains were grown on the indicated medium to measure 
the spreading of heterochromatin into the ura4'^ reporter. (£,/) ChIP analysis of H3K9me2 and Swi6 levels at the ma4'^ reporters, 
normalized to actl . The data are averages of three experiments, and error bars represent standard deviation. [F] qRT-PCR analysis of 
emc5 transcript levels, normalized to actl. Data presented are averages of three experiments, and error bars represent standard 
deviation. 



Because the effect of IRC IRA on heterochromatin 
spreading might be a result of moving the ura4^ reporter 
closer to heterochromatin, we also tested whether IRCIR 
is sufficient to block heterochromatin spreading. Ectopic 
insertion of a part of a pericentric repeat (L5) induces 
silencing of an adjacent ura4'^ reporter gene (Partridge et al. 
2002; Sadaie et al. 2008). Interestingly, when we inserted 
IRCIR between L5 and ura4'^ (Fig. IC), heterochromatin 
spreading was blocked, as indicated by the failure of cells 
to grow on FOA medium (Fig. ID). Thus, IRCIR is both 
necessary and sufficient for boundary function. 

Bdf2 is required for proper boundary function 

To identify factors required for boundary function, we 
screened a fission yeast haploid deletion library containing 
—3500 individual deletions of nonessential genes (Kim 
et al. 2010) for mutants that showed heterochromatin 



spreading (Supplemental Fig. SI). Our screen identified 
epelA, which is known to be required for boundary 
function (Ayoub et al. 2003; Zofall and Grewal 2006; 
Trewick et al. 2007), as well as an uncharacterized null 
mutant of SPAC631.02, which has been named bdf2'^ 
(Garabedian et al. 2012) and nrcl^ in PomBase. The Bdf2 
protein contains two bromo domains and is homologous to 
the mammalian bromodomain and extra terminal (BET) 
family of double bromodomain proteins (Supplemental 
Fig. S2; Belkina and Denis 2012; Prinjha et al. 2012), which 
were recently identified as therapeutic targets for a number 
of cancers (Dawson et al. 2011; Zuber et al. 2011). 

We confirmed that bdf2A resulted in heterochromatin 
spreading into IRCIR ::ura4^, as indicated by growth on 
FOA plates (Fig. IB). Furthermore, heterochromatin hall- 
marks such as H3K9 dimethylation (H3K9me2) and Swi6 
spread to the reporter gene, as indicated by chromatin 
immunoprecipitation (ChIP) analyses (Fig. IE). His tone 3 
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levels at the reporter were unaffected by hdf2A (Supple- 
mental Fig. S3 A), indicating that the change in H3K9me2 
levels Mras specific for this modification. Moreover, the 
expression level of emcS, which is the gene immediately 
adjacent to IRCIR, was reduced in bdf2A cells (Fig. IF), 
suggesting that heterochromatin spreading regulates the 
expression of neighboring genes. 

Similarly, bdf2A cells showed heterochromatin spread- 
ing into ura4'^ inserted outside of the centromere in 
heterochromatin boundary (IRC3L::ma4^] (Fig. IG-I) and 
the silent mating type region heterochromatin boundary 
iIRR::ma4^] (Supplemental Fig. S4A,B; Singh and Klar 
2008). Given that IRs and IRCs do not share sequence 
homology, our results indicate that Bdf2 is required for 
proper boundary functions at different locations, inde- 
pendent of the presence of IRCs. 

Bdf2 is enriched at heterochromatin boundaries 

To examine the cellular functions of Bdf2, we generated 
a C -terminally Flag- tagged version of Bdf2 expressed from 
its endogenous chromosomal location. Bdf2-Flag was 
functional, as no defects in boundary function were ob- 
served (Supplemental Fig. S5A). We performed ChlP-chip 
analyses of Bdf2-Flag and H3K9me2 using an Agilent 
whole-genome microarray with additional customized 
probes that provide high coverage of centromeres. Bdf2 
was enriched at all IRC elements, which correspond to 
sharp decreases in H3K9me2 levels (Fig. 2A,E). The 
localization of Bdf2 to IRCs was further confirmed by 
Chip and quantitative real-time PGR (qPGR)-based anal- 
yses (Fig. 2B,F). Bdf2 also localized at the right border of 
centromere II heterochromatin, which correlates with 
a sharp drop of H3K9me2 levels but does not contain any 
IRC element (Fig. 2G,D). Interestingly, Bdf2 was not 
enriched at the left side of the centromere II heterochro- 
matin boundary (Fig. 2C). We also found that Bdf2 local- 
ized at the mating type region boundary IR, although at 
lower levels compared with centromere boundaries (Sup- 
plemental Fig. S4G). Moreover, we found that Bdf2 also 
localized selectively at the promoter regions of a small 
group of genes (Supplemental Fig. S6A,B), indicating a 
possible role of Bdf2 in transcriptional regulation in 
addition to its role in boundary function. 

Bdf2 interacts with Epel 

To investigate further the function of Bdf2, we performed 
affinity purification of Bdf2-Flag and identified interacting 
proteins by mass spectrometry (Supplemental Fig. S7A; 
Supplemental Table SI). Significantly, among the proteins 
identified that specifically associated with Bdf2 but were 
absent in control purifications was the boundary protein 
Epel (Fig. 3 A), consistent with Bdf2's boundary function. 
Bdf2 was also associated with a number of TAFs (TBP- 
associated factors) (Supplemental Fig. S7B), which might 
explain the promoter localization of Bdf2 in euchromatic 
regions. GhIP analyses of TaflO-Flag indicated that TAFs 
are not enriched at heterochromatin boundaries (Supple- 
mental Fig. S7G). Conversely, affinity purification of Epel- 
Flag also identified Bdf2 (Fig. 3 A; Supplemental Table SI). 



To confirm the interaction between Bdf2 and Epel, we 
performed coimmunoprecipitation analysis with extracts 
from cells expressing both Bdf2-myc and Epel -Flag from 
their endogenous chromosomal loci. When Epel -Flag was 
immunoprecipitated with anti-Flag-agarose resin, Bdf2 
specifically copurified as well, demonstrating that Bdf2 
interacts with Epel in vivo (Fig. 3B). Moreover, the inter- 
action persisted when the cell lysates were treated with 
ethidium bromide (Fig. 3B), suggesting that the interaction 
between Epel and Bdf2 is direct rather than DNA-mediated. 

Epistasis analysis showed that a bdf2A epel A double 
mutant resulted in heterochromatin spreading similar to 
that of epel A, suggesting that Bdf2 and Epel function in 
the same pathway (Fig. 3C). We also consistently ob- 
served that epelA resulted in a slightly greater increase 
in heterochromatin spreading than bdf2A, suggesting 
additional roles for Epel in regulating boundary functions 
(Fig. 3C). 

To examine the relationship between Bdf2 and Epel, 
we performed ChIP analysis of Bdf2-Flag in an epelA 
strain. Interestingly, Bdf2 localization at IRCl and IRCS 
was lost or significantly reduced in epel A cells (Fig. 3D). 
Bdf2 protein levels and Bdf2 localization to the promoter 
of oca8 were not affected by epelA (Supplemental Figs. 
S6C, S8), indicating that Epel only regulates the locali- 
zation of Bdf2 at IRCs. Epel is recruited to heterochroma- 
tin by Swi6 (Zofall and Grewal 2006; Trewick et al. 2007; 
Sadaie et al. 2008; Braun et al. 201 1). As expected, Bdf2 was 
also delocalized from IRCl and IRCS in swi6A cells (Fig. 
3D). On the other hand, the localization of Epel to IRCs was 
not affected by bdf2A (Fig. 3D). These results suggest that 
Epel functions upstream of Bdf2 for localization to IRCs. 

We further demonstrated that Bdf2 and Epel interacted 
in a yeast two-hybrid assay (Fig. 3E,F), suggesting that the 
interaction between these two proteins might be direct. 
To determine the region of Bdf2 that mediates interaction 
with Epel, we generated Bdf2 fragments that lack the 
N-terminal region (AN), the C-terminal region (AC), or both 
(BD) (Fig. 3E). We found that Bdf2 lacking the C-terminal 
region failed to interact with Epel in yeast two-hybrid 
assays (Fig. 3F). Moreover, coimmunoprecipitation showed 
that Bdf2-AC-myc interaction with Epel -Flag was strongly 
reduced (Fig. 3G), suggesting that the C-terminal region of 
Bdf2 is required for Bdf2-Epel interaction in vivo. As ex- 
pected, Bdf2-AC localization to IRCl was reduced (Fig. 3H), 
concomitant with heterochromatin spreading into IRCIR:: 
uraO" (Fig. 31). 

The bromodomains of Bdf2 are required for boundary 
function 

Bromodomains are well known for their ability to bind 
acetylated lysines (Mujtaba et al. 2007). Bdf2 contains two 
bromodomains (Fig. 4A), but their binding specificities 
are unknown. Binding assays of recombinant GST-Bdf2- 
BD (amino acids 229-497, encompassing both bromodo- 
mains) with an array of his tone tail peptides containing 
different combinations of post-translational modifica- 
tions showed that Bdf2 preferentially binds to multiply 
acetylated histone H4 tail peptides (Supplemental Fig. S9). 
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Figure 2. Bdf2 is enriched at IRCs. {A,C,E) ChlP-chip analysis of Bdf2-Flag and H3K9me2 levels across centromeres I, II, and III. The 
data are averages of two microarrays. Black arrows in the diagrams indicate the positions of IRCs. Bars indicate tRNA genes. Note that 
centromere II does not contain IRCs. {B,D,F] ChlP-qPCR analysis of Bdf2-Flag levels at boundary regions, normalized to actl. The 
positions of amplified regions used for qPCR quantification are indicated. The data are averages of three experiments, and error bars 
represent standard deviation. 



Such results are consistent v^ith those of a systematic 
study of BET family bromodomains, which indicates that 
they usually have broad specificity and preferentially bind 
to multiply acetylated his tone tails (Filippakopoulos et al. 
2012). To confirm this finding, we generated recombinant 
Fiis-Bdf2-BD and examined its binding to a tetra-acetylated 
histone Fi4 tail peptide (K5, K8, K12, and K16) and a 
diacetylated histone H3 peptide (K9 and K14). Bdf2-BD 
bound strongly to the acetylated Fi4 peptide but not to its 
unacetylated form (Fig. 4B). Moreover, mutations of two 
amino acids that are predicted to form acetyl-lysine-binding 
sites, one in each of the bromodomains (Y268A and Y430A, 



hereafter denoted as 2YA) (Supplemental Fig. S2; Mujtaba 
et al. 2007), abolished this binding, indicating that Bdf2 
binds to acetylated histone tails through its bromodomains 
(Fig. 4B). Bdf2 had little affinity for the diacetylated 
histone Fi3 peptide or any of the singly acetylated Fi4 
peptides (Fig. 4B; data not shown), so it seems to prefer 
multiply acetylated Fi4 tails. 

To examine further the role of bromodomains of Bdf2 in 
boundary function, we introduced the 2YA mutations at 
the endogenous bdf2'^ chromosomal location. Interest- 
ingly, boundary function was compromised to a degree 
similar to that of bdf2A, as indicated by silencing of 
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Figure 3. Epel recruits Bdf2 to IRCs. [A] Multidimensional protein identification technology (MudPIT) mass spectrometry analyses of 
purified Bdf2-Flag and Epel -Flag complexes. The number of Bdf2 and Epel peptides identified and the percentage of each protein that 
these peptides cover are indicated. {B,G) Cell extracts from the indicated strains were incubated with Flag-agarose beads to isolate Epel- 
Flag. Bound proteins were resolved by SDS-PAGE, and Western blot analyses were performed with Myc and Flag antibodies. (C,7) Serial 
dilution analyses were performed to measure heterochromatin spreading outside oilRClR:: ura4'^. {D,H] ChIP analysis of Bdf2-Flag and 
Epel -Flag levels at IRCs, normalized to actl. The data are averages of three experiments, and error bars represent standard deviation. {E] 
Diagram of Bdf2 constructs used in yeast two-hybrid analysis. [F] Yeast two-hybrid analysis of Bdf2 with Epel. Bdf2 was fused with the 
Gal4 DNA-binding domain, and Epel was fused with an activation domain. Interaction between Bdf2 and Epel resulted in the 
activation of a HISS reporter gene, allowing cells to grow in the absence of histidine. 



IRCIR ::um4^ (Fig. 4C) and the spreading of H3K9me2 and 
Swi6 into the reporter gene (Fig. 4D). These results suggest 
that binding of acetylated histone Fi4 tails by Bdf2 is 
critical for boundary function. The mutations had no effect 
on Bdf2 protein levels (Fig. 4F), and Bdf2 still localized at 
IRCs, although at lower levels compared with that of wild 
type (Fig. 4E). This result suggests that binding of Bdf2 to 
acetylated histones also stabilizes Bdf2 at IRCs. Neither 
the Y268A nor the Y430A single mutant affected hetero- 



chromatin boundary function (Supplemental Fig. S5B), 
suggesting that the two bromodomains have redundant 
roles, at least in regulating boundary formation at IRCs. 

Tethering Bdf2 to DNA is sufficient to establish 
a heterochromatin boundary 

We further tested whether recruiting Bdf2 was sufficient 
to establish a heterochromatin boundary independent of 
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Figure 4. The bromodomains of Bdf2 are required for boundary function. [A] Schematic diagram of Bdf2 protein. The positions of 
mutations are indicated. [B] Peptide pull-down assays of the indicated recombinant proteins with biotinylated H4 tail (2-24) unmodified 
and tetra-acetylated at the K5, K8, K12, and K16 (Ac-H4) peptides as well as H3 tail (1-21) unmodified and diacetylated at the K9 and 
K14 (Ac-H3) peptides. (C) Serial dilution analysis to measure heterochromatin spreading outside of IRClR::ura4^. {D,E] ChIP analysis 
of H3K9me2 and Swi6 levels at IRClR::ura4^ and Bdf2-Flag levels at IRCIR and IRC3L, normalized to actl. The data are averages of 
three experiments, and error bars represent standard deviation. [F] Western blot analysis of Bdf2-Flag levels. (G) Schematic diagram of 
the L5-gbs-uia4^ reporter. (H) Serial dilution analysis to measure heterochromatin spreading to ma4'^. All strains used contain epelA to 
rule out the possibility that Bdf2 recruits Epel to establish heterochromatin boundaries. (7) ChIP analysis of H3K9me2 and Swi6 levels 
at L5-gbs-ura4^ and Bdf2-GBD levels at 3xgbs, normalized to actl. The data are averages of three experiments, and error bars represent 
standard deviation. All strains used contain epelA. 



Epel. We inserted three copies of the Gal4-binding site 
(gbs) between L5 and ma4'^ to create the L5-gbs-ma4'^ 
reporter (Fig. 4G). In the absence of Epel, Bdf2-GBD, but 
not Bdf2-2YA-GBD, was able to block the spreading of 
heterochromatin from L5 to ura4'^, as indicated by silenc- 
ing assays (Fig. 4Fi) and ChIP analyses of Fi3K9me2 and 
Swi6 (Fig. 41). Bdf2-GBD and Bdf2-2YA-GBD were equally 
enriched at gbs (Fig. 41), indicating that it is not the 
binding of Bdf2 per se but rather its ability to bind 
acetylated histone Fi4 tails that is essential for boundary 
activity. 



Bdf2 counteracts Sir 2 -mediated histone deacetylation 

The spreading of heterochromatin from nucleation sites 
depends on not only Clr4-mediated Fi3K9me and binding 
of Swi6, but also Sir2-mediated histone deacetylation 
(Fiall et al. 2002; Shankaranarayana et al. 2003; Zhang 
et al. 2008; Buscaino et al. 2013;). Fiowever, sir2A has only 
minor effects on silencing of reporter genes inserted at 
pericentric heterochromatin (Shankaranarayana et al. 2003; 
Freeman-Cook et al. 2005; Alper et al. 2013; Buscaino et al. 
2013), and it is not known whether Sir2 is required for 



1892 GENES & DEVELOPMENT 



Heterochromatin boundary formation by Bdf2 



heterochromatin spreading outside IRCs. We generated 
bdf2A sii2A and epelA sir2A strains with the IRCIR:: 
VLia4^ reporter. In both cases, heterochromatin failed to 
spread to the reporter gene, as indicated by silencing assays 
and Chip analyses (Fig. 5A-D), suggesting that heterochro- 
matin spreading outside of its boundaries depends on Sir2. 

The role of H4K16ac in regulating heterochromatin 
spreading 

The major substrates of Sir2 are H3K9ac and H4K16ac in 
vitro (Shankaranarayana et al. 2003; Alper et al. 2013). 
Since Clr4-mediated H3K9me is critical for heterochro- 
matin assembly (Rea et al. 2000; Nakayama et al. 2001), it 
is conceivable that Sir2 functions to generate deacety- 
lated H3K9 for subsequent action by Clr4. However, 
whether deacetylation of H4K16 is also involved in hetero- 
chromatin spreading is unknown. The fact that Bdf2 and 
Epel counteract Sir2 for boundary function and that Bdf2 
binds to acetylated histone H4 instead of H3 tail suggests 



that H4K16ac also plays a role in regulating heterochro- 
matin spreading. Indeed, ChIP analysis showed that 
H4K16ac levels at IRCl were significantly reduced in 
epel A and bdf2A cells (Fig. 5E). Moreover, H4K16ac levels 
were increased in sir2A cells but restored to near wild-type 
levels in epel A sir2A and hdf2A sir2A cells (Fig. 5E). 
Furthermore, recombinant Sir2 showed robust activity in 
an in vitro histone deacetylase (HDAC) assay when a tetra- 
acetylated histone H4 peptide was used as substrate, and 
this activity was strongly reduced when recombinant Bdf2, 
but not Bdf2-2YA, was added to the reaction (Fig. 5F). 

If the function of Bdf2 is to protect H4K16ac, we expect 
that H4K16ac levels will be higher at heterochromatin 
boundaries where Bdf2 is present. To test this idea, we per- 
formed ChlP-chip analyses of H4K16ac at heterochromatin 
regions. H4K16ac levels were low in the middle of heterochro- 
matin, which is consistent with the fact that heterochro- 
matin is generally devoid of histone acetylation. Interest- 
ingly, however, H4K16ac levels spiked at heterochromatin 
boundaries, corresponding to Bdf2 peaks (Fig. 6A-C). 
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IRC1R::ura4^ control FOA IRC1R::ura4^ control FOA 




Figure 5. Bdf2 counteracts Sir2-mediated histone deacetylation. (A,C) Serial dilution analysis to measure heterochromatin spreading 
outside of IRCIR •.•.ma4^. {B,D,E) ChIP analysis of H3K9me2 and Swi6 levels at IRCIR ::ma4^ and H4K16ac levels at IRCl, normalized 
to actl. The data are averages of three experiments, and error bars represent standard deviation. {F, left] Coomassie-stained gel of 
recombinant Sir2. [Right] HDAC assays were performed with recombinant Sir2, a tetra-acetylated histone H4 peptide, and recombinant 
Bdf2. The production of nicotinamide was measured via the generation of nicotinic acid and free ammonia by nicotinamidase. 
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Figure 6. H4K16ac modulates heterochromatin boundary function. [A-C] ChlP-chip analysis of H4K16ac levels across centromeres I, 
II, and III. The Bdf2-Flag profile is also shown for comparison. [D] Serial dilution analysis to measure heterochromatin spreading outside 
oilRClR:: ura4^. All strains used are in an h4. lA/h4.3A background. {E,F) ChIP analysis of H3K9me2 and Swi6 levels at IRCIR ::ura4^ 
and Bdf2 levels at IRCl, normalized to actl. The data are averages of three experiments, and error bars represent standard deviation. All 
strains used are in an h4.1A/h4.3A background. 



To test directly whether deacetylation of H4K16 is re- 
quired for heterochromatin spreading in fission yeast, we 
generated an H4K16Q mutation, which mimics H4K16ac. 
Fission yeast contains three pairs of H3/H4 genes, and 
deleting two of the Fi3/Fi4 pairs {h4.1A/h4.3A) allows the 
production of histone Fi4 from the remaining single copy 
(Mellone et al. 2003). In this background, histone Fi4 
levels were reduced, as indicated by Western blot analyses 
(Supplemental Fig. SI OA). Fie tero chromatin still spread 
into IRCIR ::ura4'^ in epelA cells albeit at lower efficiency. 



indicating that histone dosage also contributes to hetero- 
chromatin spreading (Supplemental Fig. SI OB). Interest- 
ingly, in bdf2-2YA H4K16Q cells, heterochromatin failed 
to spread into IRCIR ::ura4^, as indicated by serial di- 
lution analyses and ChIP analyses of FF3K9me2 and Swi6 
(Fig. 6D,E), suggesting that deacetylation of FF4K16 is in- 
deed required for heterochromatin spreading. Moreover, a 
H4K16R mutation, which abolished FF4K16ac, resulted in 
heterochromatin spreading into IRCIR ::ma4'^ (Fig. 6D,E). 
Similar results were obtained with epelA and bdf2A 
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(Supplemental Fig. SI 1). The localization of Bdf2 to IRCl 
was reduced in H4K16R cells but slightly increased in 
H4K16Q cells (Fig. 6F). In contrast, the localization of 
Bdf2-2YA to IRCl was little affected by H4K16 mutations 
(Fig. 6F). These results demonstrate that Bdf2 is critical 
for regulating H4K16ac levels at boundaries and hetero- 
chromatin spreading. 

Mstl regulates H4K16ac levels at IRCs 

Among the histone acetyltransf erases (HATs) in fission 
yeast, Mstl is a member of the highly conserved MYST 
family of HATs, and it forms a complex similar to the 
NuA4 complex of budding yeast (Gomez et al. 2005, 2008; 
Shevchenko et al. 2008), which acetylates H4K16 (Lee and 
Workman 2007). We found that a mstl-1 allele, which 
contains a temperature-sensitive L344S mutation (Gomez 
et al. 2008), resulted in heterochromatin spreading past 
IRCIR even at a permissive temperature, as indicated by 
the silencing of IRCIR ::ura4^ and the accumulation of 
H3K9me2 and Swi6 at the reporter (Fig. 7A,B). Moreover, 
in mstl-1 H4K16Q cells, heterochromatin spreading is 
abolished (Fig. 7A,B). Furthermore, mstl-1 resulted in a 
significant decrease of H4K16ac levels and reduced re- 
cruitment of Bdf2 at IRCl (Fig. 7C,D). Thus, Mstl directly 
regulates H4K16ac levels and the recruitment of Bdf2 at 
IRCs for proper boundary formation. 

Discussion 

Proper formation of boundaries between chromosomal 
domains is essential for maintaining stable gene expres- 
sion patterns (Gaszner and Felsenfeld 2006; Valenzuela 
and Kamakaka 2006). Here we showed that a double 
bromodomain protein of the BET family, Bdf2, is both 
necessary and sufficient for the formation of boundaries 
that limit the spreading of heterochromatin. We further 
showed that the bromodomains of Bdf2, which preferen- 
tially bind to multiply acetylated histone H4 tails, are 
critical for Bdf2 boundary function. Moreover, we dem- 
onstrated that Bdf2 directly protects histones from Sir2- 
mediated deacetylation and that the deacetylation of 
H4K16 is necessary for H3K9me and Swi6-mediated het- 
erochromatin spreading. Below, we discuss the implications 
of our findings for heterochromatin spreading and boundary 
functions. 

The role of Bdf2 in boundary formation 

The formation of heterochromatin has long been regarded 
as a paradigm for the assembly of self-propagating chro- 
matin structures (Moazed 2011). The mechanism of 
heterochromatin spreading is complex and not well un- 
derstood (Talbert and Henikoff 2006). A simplified view is 
that heterochromatin spreads by a positive feedback loop 
of histone modification and chromatin protein binding, 
the reiteration of which leads to the spreading of these 
modifications/proteins from nucleation sites (Talbert and 
Henikoff 2006; Grewal and Elgin 2007). There are two 
ways heterochromatin spreading is curbed. One is through 
competition between positive and negative forces, which 
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Figure 7. {A] Serial dilution analysis to measure heterochro- 
matin spreading outside of IRCIR ::ura4'^. All strains used are in 
an h4.1A/h4.3A background. [B-D] ChIP analysis of H3K9me2 
and Swi6 levels at IRCIR ::ma4^ and H4K16ac and Bdf2 levels 
at IRCl, normalized to actl. The data are averages of three 
experiments, and error bars represent standard deviation. (£) A 
model by which H4K16ac and Bdf2 regulate heterochromatin 
spreading and boundary function. Heterochromatin spreading is 
mediated by cycles of Clr4-mediated H3K9me and the recruit- 
ment of Swi6. Sir2 mediates the deacetylation of H3K9, which is 
required for creating the substrate for Ch4, as well as the 
deacetylation of H4K16, which facilitates nucleosome compac- 
tion to bring the adjacent nucleosome closer to Clr4. Bdf2 is 
recruited to IRCs through Epel, which is highly enriched at 
IRCs. Mstl mediates H4K16ac at IRCs, which further stabilizes 
binding of Bdf2 to IRCs through its bromodomains. Bdf2 pro- 
tects H4 tails from deacetylation by Sir2, preventing further 
heterochromatin spreading. 



GENES & DEVELOPMENT 1895 



Wang et al. 



results in the formation of extended transition zones that 
fluctuate when the heterochromatin-euchromatin bal- 
ance changes (Kimura and Horikoshi 2004). These borders 
are termed ^'negotiable borders/' since they are not defined 
by DNA sequence but move depending on the dosage of 
opposing activities. For example, in the classical example 
of position effect variegation (PEV) in Drosophila, where 
the white gene is juxtaposed to pericentric heterochroma- 
tin due to chromosome rearrangement, it is stochastically 
silenced by heterochromatin in a portion of cells due to 
variegation in the distance of heterochromatin spreading, 
resulting in mottled eyes (Elgin and Renter 2007). Most 
of the genes identified that affect PEV are involved in 
regulating chromatin function in general rather than 
boundary functions specifically (Elgin and Renter 2007). 
The more reliable way to stop heterochromatin from spread- 
ing is by specific DNA elements that actively establish 
borders to stop the spreading of heterochromatin (Gaszner 
and Felsenfeld 2006; Valenzuela and Kamakaka 2006). 
These DNA elements recruit histone-modifying enzymes 
to directly modify histones to counteract heterochroma- 
tin spreading (Old et al. 2004; West et al. 2004). Alterna- 
tively, they tether chromosomal regions to the nuclear 
periphery, resulting in physical separation of chromosomal 
domains (Ishii et al. 2002; Noma et al. 2006). Regardless of 
mechanism, the key to boundary formation is to break up 
the self-propagation of heterochromatic histone modifi- 
cations. We showed here that Bdf2 disrupts this chain 
reaction by protecting modified histone tails from hetero- 
chromatic modifications through the recruitment of a his- 
tone tail-binding protein to the boundary region (Fig. 7E). 

Bdf2 is recruited to IRCs through boundary protein Epel 
(Fig. 3D). This recruitment is dependent on the C-terminal 
region of Bdf2, which mediates its interaction with Epel 
(Fig. 3F-H). Loss of Epel resulted in the delocalization of 
Bdf2 from IRCs without affecting the localization of Bdf2 
to gene promoters elsewhere in the genome (Fig. 3D; 
Supplemental Fig. S6C). However, Epel exhibits a broad 
localization pattern across entire heterochromatin do- 
mains, although Epel levels at IRC elements are higher 
compared with those in the middle of heterochromatin 
(Zofall and Grewal 2006; Braun et al. 201 1). This suggests 
that high concentrations of Epel are required for Bdf2 
localization to IRCs. Alternatively, other factors present 
at IRCs might enhance the interaction between Epel and 
Bdf2 to promote the localization of Bdf2 to IRCs. We favor 
the latter idea because even in ddhlA cells, which exhibit 
a significant increase of Epel at the body of heterochro- 
matin, there is only marginal enrichment of Bdf2 at the 
heterochromatic repeats (Supplemental Fig. SI 2). It is 
possible that acetylated histone H4 stabilizes Bdf2 bind- 
ing to chromatin, as indicated by lower levels of Bdf2 at 
IRCs when the bromodomains are compromised or when 
H4K16ac levels are reduced (Figs. 4E, 6F, 7D). The lack of 
Bdf2 enrichment at pericentric repeats in ddhlA cells 
is also consistent with this idea, since histones within 
heterochromatin regions are generally hypoacetylated 
(Kurdistani and Grunstein 2003; Grewal and Elgin 2007). 
The localization of Bdf2 at the right side of the centromere 
II heterochromatin boundary (Fig. 2B,E), which does not 



show elevated Epel levels (Zofall and Grewal 2006), is 
also consistent with the idea that additional mechanisms 
contribute to Bdf2 recruitment to heterochromatin bound- 
aries. The high concentration of Epel at heterochromatin 
boundaries is achieved through Cul4-Ddbl -mediated 
degradation of Epel in the middle of heterochromatin 
(Braun et al. 2011) and thus should be independent of 
DNA sequences. However, heterochromatin borders at 
centromeres and the silencing mating type region are 
well defined, suggesting that boundary DNA elements 
have independent contributions to shape the localization 
of Epel and Bdf2. 

Overexpression of Epel results in defects in pericentric 
heterochromatin assembly (Zofall and Grewal 2006; 
Trewick et al. 2007) even in hdf2A cells (Supplemental Fig. 
S13), suggesting that Epel negatively regulates hetero- 
chromatin assembly independently of its role in recruiting 
Bdf2. Such an activity may also contribute to the boundary 
function of Epel. This idea is further supported by the fact 
that epelA showed stronger spreading of heterochromatin 
than hdf2A. Although Epel has also been proposed to be 
a demethylase or a hydroxylase and mutations of putative 
active site residues affect Epel function (Trewick et al. 
2005, 2007), the substrate of Epel has not been identified. 
Further characterization of the enzymatic activity of Epel 
is needed to reveal additional mechanisms by which Epel 
regulates heterochromatin spreading. 

In budding yeast, the double bromodomain protein 
Bdfl is also required to prevent Sir2-mediated hetero- 
chromatin spreading near telomeres through the binding 
of Bdfl to acetylated histone tails (Ladurner et al. 2003). 
Unlike fission yeast Bdf2, budding yeast Bdfl does not 
specifically localize to defined boundary elements but is 
recruited to the transition zone, in which a gradient of 
H4K16ac is formed by the actions of HAT SAS2 and 
HDAC Sir2 (Kimura et al. 2002; Suka et al. 2002; Ladurner 
et al. 2003). Bdfl is also part of the Swrl complex that is 
required for the deposition of histone variant H2A.Z 
(Krogan et al. 2003; Kobor et al. 2004; Mizuguchi et al. 
2004), which is also involved in counteracting heterochro- 
matin spreading (Meneghini et al. 2003). There is a Bdfl 
homolog in fission yeast, which is part of the fission yeast 
Swrl complex (Buchanan et al. 2009; Hou et al. 2010). 
However, fission yeast Bdfl is not involved in boundary 
function at IRCs (Supplemental Fig. S14). 

H4K16ac and heterochromatin spreading in fission 
yeast 

The bromodomains of Bdf2 interact with tetra-acetylated 
but not singly acetylated histone H4 tail peptides in vitro. 
This is consistent with studies of human BET family 
proteins, which favor multiply acetylated substrates 
(Filippakopoulos et al. 2012). Mutational analysis of acet- 
ylated lysines on the histone H4 tail in budding yeast 
demonstrates that Lys5, Lys8, Lysl2, and Lysl6 are par- 
tially redundant and cumulatively regulate gene expres- 
sion, heterochromatin assembly, maintenance of genome 
integrity, and cell cycle progression, although Lysl6 acet- 
ylation has a more dominant role (Megee et al. 1990, 1995; 
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Park and Szostak 1990; Durrin et al. 1991; Dion et al. 
2005). The fact that fission yeast Sir2 preferentially deacety- 
lates K16 among the H4 tail lysines (Shankaranarayana 
et al. 2003; Alper et al. 2013) and that Bdf2 counteracts Sir2 
to regulate H4K16ac levels at IRCs suggests that acety- 
lated H4K16 is a major target of Bdf2 in vivo. Indeed, the 
association of Bdf2 with IRCs and gene promoters is 
significantly reduced in H4K16R cells (Fig. 6F; Supple- 
mental Fig. S6D). 

H4K16ac directly regulates the compaction of nucleo- 
somal arrays in vitro (Shogren-Knaak et al. 2006), and in 
budding yeast, H4K16ac is critical for regulating hetero- 
chromatin spreading in vivo (Kimura et al. 2002; Suka 
et al. 2002). In fission yeast, eliminating H4K16ac with 
a H4K16R mutation resulted in heterochromatin spread- 
ing outside of IRCIR, and mimicking hyperacetylation of 
H4K16 with a H4K16Q mutation effectively blocked 
heterochromatin spreading in bdf2A and epelA cells. 
Such results suggest that deacetylation of H4K16 is an 
integral part of H3K9me and Swi6-mediated heterochro- 
matin spreading. However, H4K16ac-mimicking muta- 
tions such as H4K16A and H4K16Q have little effect on 
silencing of reporter genes inserted inside pericentric 
repeats such as otr::ade6^ (Supplemental Fig. SI 5; Mellone 
et al. 2003), indicating that additional heterochromatin 
assembly pathways at pericentric regions, such as RNAi, 
can overcome the requirement of Sir2 (Alper et al. 2013; 
Buscaino et al. 2013). 

In both budding yeast and fission yeast, Sir2 is the 
conserved HDAC that regulates heterochromatin assem- 
bly (Kimura et al. 2002; Suka et al. 2002; Shankaranarayana 
et al. 2003; Freeman-Cook et al. 2005; Buscaino et al. 2013). 
In fission yeast, it was proposed that Sir2-mediated 
deacetylation of H3K9 is a prerequisite for H3K9me 
(Shankaranarayana et al. 2003). However, sii2A has only 
minor effects on H3K9me levels at heterochromatin 
nucleation centers such as pericentric repeats or the cenH 
sequence at the mating type region (Shankaranarayana 
et al. 2003; Freeman-Cook et al. 2005; Alper et al. 2013; 
Buscaino et al. 2013), suggesting that deacetylation of 
H3K9 in the absence of Sir2 can be accomplished by other 
HDACs. Our results that Sir2 regulates H4K16ac levels at 
heterochromatin boundaries and that H4K16 mutants 
modulate heterochromatin boundary function suggest 
that Sir2 also deacetylates H4K16 in regulating hetero- 
chromatin spreading. We hypothesize that the deacetyla- 
tion of H4K16 increases chromatin compaction levels, 
thus bringing the adjacent nucleosomes closer to Clr4 for 
repeated cycles of H3K9me propagation. The fact that 
reducing histone dosage affects the ability of heterochro- 
matin to spread (Supplemental Fig. SIO) also supports 
such a hypothesis. 

In budding yeast, Sas2 is the major H4K16 acetyltrans- 
ferase that regulates silencing and boundary formation 
(Kimura et al. 2002; Suka et al. 2002; Shia et al. 2006). 
However, there is no Sas2 homolog in fission yeast. We 
discovered that the fission yeast MYST family acetyl- 
transferase Mstl is required for H4K16ac and the binding 
of Bdf2 at IRCs, and a mutation in Mstl also resulted in 
defects in boundary formation at IRCs. Mstl forms a 



complex with composition similar to that of the budding 
yeast NuA4, which acetylates H4K16 as well as other 
residues (Lee and Workman 2007; Gomez et al. 2008; 
Shevchenko et al. 2008). However, in budding yeast, 
mutations in Esal, the catalytic subunit of NuA4, have 
no effect on telomeric heterochromatin assembly (Kimura 
et al. 2002; Suka et al. 2002). Thus, the fission yeast Mstl 
complex has taken over the role of both Sas2 and NuA4 
complexes in budding yeast, which might be attributed to 
the slight difference in the composition of Mstl and 
NuA4 complexes (Shevchenko et al. 2008). 

Ovei expression of Swi6 enhances heterochromatin 
spreading 

Multiple mechanisms cooperate to establish heterochro- 
matin boundaries. In addition to Epel, TFIIIC binding to 
IR elements contributes to boundary formation at the 
silent mating type region, and tRNA genes and transcrip- 
tion oilRCl contribute to boundary activity at centromeres 
(Noma et al. 2006; Scott et al. 2006; Keller et al. 2013). As 
a result, heterochromatin spreading in epel A cells is a 
stochastic event happening only in a small portion of 
cells. Although the reporter-based assays provide sensi- 
tive readouts of heterochromatin spreading in this sub- 
population, it is difficult to perform biochemical assays 
with ensembles of cells, since the majority of cells do not 
show spreading. For example, our ChlP-chip analysis 
showed that in epelA cells without Swi6 overexpression 
or reporter-based selection, H3K9me2 levels around Ji^Ci 
were similar to those of wild- type cells (data not shown), 
consistent with earlier analysis of H3K9me2 levels in epelA 
cells at the mating type region and centromeres (Ayoub 
et al. 2003; Braun et al. 201 1). The use of an additional copy 
of Swi6 provides a more profound effect on heterochromatin 
spreading at centromere boundaries, and similar approaches 
have been used to enhance heterochromatin spreading at 
the mating type region (Noma et al. 2001, 2006). 

The functions of Bdf2 in other processes 

In addition to its role in boundary function, Bdf2 also 
interacted with a number of TAFs, and our ChlP-chip 
analysis showed that Bdf2 localized at a small group of 
gene promoters (—13% of all promoter probes in our 
microarray), suggesting that Bdf2 might be involved in 
transcriptional regulation of certain genes. IRCs were 
transcribed (Zofall and Grewal 2006; Keller et al. 2013), 
and IRCl transcript levels were strongly reduced in bdf2A 
and epel A cells (Supplemental Fig. SI 6), suggesting that 
Bdf2 regulates transcription of IRCs. However, targeting 
of Bdf2 to Gal4-binding sites is sufficient to block the 
spreading of heterochromatin at an ectopic site in a bro- 
modomain-dependent manner, suggesting that the ability 
of Bdf2 to counteract heterochromatin spreading is 
through protection of H4K16ac by its bromodomains, 
although transcription of IRCs has an independent con- 
tribution to boundary function (Keller et al. 2013). It is 
also possible that Bdf2 regulates transcription of other 
factors involved in heterochromatin boundary functions, 
and it would be interesting to determine the role of Bdf2 
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in regulating gene expression at euchromatic regions in the 
future. Bdf2 has also been shown to regulate DNA damage 
response, which is attributed to its effect on regulating 
global H4 acetylation levels (Garabedian et al. 2012). 

Implications of BET protein functions in other systems 

The BET family of bromodomains in mammals includes 
BRD2, BRD3, BRD4, and BRDT. They associate with 
diverse protein complexes involved in chromatin modi- 
fications and play important roles in transcription regu- 
lation, cell cycle control, gene bookmarking, and viral 
genome transcription (Belkina and Denis 2012). Misregu- 
lation of BET proteins, especially Brd4, has been linked to 
a number of human cancers. For example, chromosome 
translocation between BRD4 and NUT (nuclear protein 
in testis) leads to a highly lethal form of carcinoma, and 
BRD4 is a therapeutic target in acute myeloid leukemia 
(French et al. 2003; Zuber et al. 2011). The identification 
of small molecules such as JQl and IBET that inhibit 
the binding of BET family bromodomains to chromatin 
(Filippakopoulos et al. 2010; Nicodeme et al. 2010) makes 
them the focus of drug-mediated epigenetic treatment of 
diverse types of cancers (Dawson et al. 2011; Delmore 
et al. 2011; Mertz et al. 2011; Zuber et al. 2011; Loven 
et al. 2013). Although the precise mechanism of these 
inhibitors is not well understood, the displacement of BET 
proteins from chromatin in part down-regulates oncogene 
c-Myc to promote differentiation (Dawson et al. 2011; 
Mertz et al. 2011; Zuber et al. 2011; Loven et al. 2013). 
However, mice treated with these drugs do not exhibit 
tissue-renewing defects associated with lower Myc activ- 
ity, indicative of additional mechanisms (Filippakopoulos 
et al. 2010; Dawson et al. 2011; Delmore et al. 2011; 
Mertz et al. 2011; Zuber et al. 2011). The fact that BET 
domain proteins Bdf2 in fission yeast and Bdf 1 in budding 
yeast regulate heterochromatin spreading suggests that 
BET inhibitors may affect the spreading of heterochro- 
matin and reset the epigenetic landscape in cancer cells. 
In other multicellular organisms such as zebrafish, worms, 
and flies, BET proteins are essential for the determination 
of cell fate during development (Huang and Dawid 1990; 
Chang et al. 2007; Dibenedetto et al. 2008; Shibata et al. 
2010). In particular, the Drosophila BET protein FSH 
(female sterile 1 homeotic) is classified as a Trithorax group 
gene that counteracts Polycomb group protein-mediated 
gene silencing (Chang et al. 2007; Kockmann et al. 2013). 
Thus, it seems that the function of BET proteins in 
antagonizing gene silencing is highly conserved across 
species. 

Materials and methods 

Fission yeast strains and genetic analyses 

The IRClR::ura4'^ reporter was constructed by inserting ura4^ 
at the right side of IRCIR (chromosome I, 3790595). The query 
strain used for our screen contains a closely linked natMX6 
cassette (chromosome I, 3804996). IRCIRA removed chromo- 
some I: 3789581-3790711. The IRC3L::ma4^ reporter was con- 
structed by inserting ura4'^ at the left side of IRC3L (chromosome 



m, 1067955). L5-IRClR-ura4^ and L5-3gbs-ura4'^ strains were 
constructed by inserting IRCIR or Sgbs sequences between L5 
and ura4'^ using L5-ura4'^ (Sadaie et al. 2004) as a template. Bdf2- 
Flag, Bdf2-myc, and Epel-Flag were constructed by a PCR-based 
module method. bdf2A, epelA, and sii2A were derived from the 
Bioneer fission yeast deletion library, verified via PGR, and back- 
crossed. The bdf2-A C and bdf2-2YA strains were constructed by 
integrating a PGR fragment containing the mutations and a Flag- 
kanMX6 cassette into the endogenous bdf 2'^ locus. H4K16R and 
H4K16Q mutations were constructed and introduced into the 
li4.1A/li4.3A background as described previously (Mellone et al. 
2003), with the exception of introducing a linked natMX6 cassette. 
Genetic crosses were used to construct all other strains. For serial 
dilution plating assays, 10-fold dilutions of a log-phase culture were 
plated on the indicated medium and grown for 3 d at 30°G. All 
strains used in heterochromatin spreading assays contain an extra 
copy of the swi6'^ gene inserted at the ade6 locus and driven by the 
ade6 promoter to enhance heterochromatin spreading. A list of 
the yeast strains used is provided in Supplemental Table S3. 

Chip analysis 

GhIP analyses were performed as described previously (Hou et al. 
2010). Antibodies used were H3K9me2 (Abeam, abl220), H4K16ac 
(Active Motif, 39167), H3 (Abeam, abl791), H4 (Abeam, abl0158), 
and Flag (Sigma, A2220). Swi6 antibody was custom-made with 
full-length recombinant Swi6 protein and affinity-purified (Reddy 
et al. 2011). 

GhlP-chip analysis was performed according to the Agilent 
Yeast GhlP-on-chip Analysis protocol. The microarray used was 
an Agilent Schizosaccharomyces pombe Whole-Genome GhlP- 
on-chip microarray (G4810A) with additional probes that en- 
compass centromeres, which were originally absent from the 
array due to the repetitive nature of these DNA sequences. 
Blunt-end DNA was generated from immunoprecipitated chro- 
matin fractions (GhIP) or whole-cell extract (WGE) with T4 DNA 
polymerase and then ligated to a linker. GhIP and WGE DNA 
were amplified from blunt-end DNA samples with primers 
annealing to the linker and were labeled by Gy5-dUTP or Gy3- 
dUTP, respectively with random priming PGR (Invitrogen GGH 
kit). Three micrograms of Gy 5 -labeled GhIP DNA and the 
corresponding Gy3-labled WGE DNA were hybridized to the 
microarray. The slides were washed and processed in accordance 
with Agilent protocols and scanned with an Agilent scanner. Data 
were collected with the Agilent Feature Extraction program. The 
enrichment value for each probe was calculated by dividing 
normalized GhIP signal by WGE signal. The data presented are 
averages of two independent GhlP-chip experiments. GhlP-chip 
data have been submitted to the Gene Expression Omnibus 
under accession number GSE46430. 

qPGR was performed with Maxima SYBR Green qPGR Master 
Mix (Thermo Scientific) in a StepOne Plus Real-Time PGR 
System (Applied Biosystems). DNA serial dilutions were used 
as templates to generate a standard curve of amplification for 
each pair of primers, and the relative concentrations of target 
sequence and a control actl sequence in GhIP and WGE samples 
were calculated accordingly. The final enrichment was calcu- 
lated as [(GhIP target)/(WGE target)]/[(GhIP actl]/{WCE actl]]. A 
list of the primers used is provided in Supplemental Table S4. 

Protein purification, coimmunoprecipitation, and mass 
spectrometry analysis 

Exponentially growing yeast cells were harvested, washed with 
2X HG buffer (300 mM HEPES-KOH at pH 7.6, 2 mM EDTA, 
100 mM KGl, 20% glycerol, 2 mM DTT, protease inhibitor 
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cocktail [Roche]), and frozen in liquid nitrogen. Crude cell ex- 
tracts were prepared by vigorously blending frozen yeast cells 
with dry ice using a household blender, followed by incubation 
with 30 mL of 1 X HC buffer containing 250 mM KCl for 30 min. 
The lysate was cleared by centrifugation at 82,700g for 3 h. The 
supernatants were precleared with protein A agarose, then 
incubated with 200 jxL of Flag-agarose overnight, and washed 
eight times with Ix HC containing 250 mM KCl. For mass 
spectrometry analysis, bound proteins were eluted with 200 jxg/mL 
3xFlag peptides followed by TCA precipitation. Multidimensional 
protein identification technology (MudPIT) mass spectrometry 
analysis was performed as described previously (Wang et al. 
2009). For coimmunoprecipitation analysis, bound proteins were 
resolved by SDS-PAGE followed by Western blot analyses with 
Myc (Santa Cruz Biotechnology, sc-789) and Flag (Sigma, F7425) 
antibodies. 

Peptide-binding assays 

The bromodomains of Bdf2 (amino acids 229-497) were cloned 
into pGEX or pRSET bacterial expression vector. The 2YA 
mutant was generated by site-directed mutagenesis according 
to the manufacturer's protocol (Agilent). Recombinant GST and 
His-tagged Bdf2 proteins were purified with Glutathione Sepharose 
4B (GE) and Talon Metal Affinity Resin (Clontech), respec- 
tively, according to manufacturer's protocols, followed by gel 
filtration with a Superdex 200 column. Binding of GST-Bdf2-BD 
to a Modified Histone Peptide Array (Active Motif, catalog no. 
17-0756-01) was performed according to the manufacturer's pro- 
tocol. The array contains 59 different post-translational modifi- 
cations for histone acetylation, methylation, phosphorylation, 
and citruUination on the N-terminal tails of histones H2A, H2B, 
H3, and H4 (Supplemental Table S2). Each 19mer peptide may 
contain up to four modifications. Array Analyze software (Active 
Motif) was used to analyze spot intensity from the image. The 
results were graphed as a specificity factor, which is the ratio of 
the average intensity of all spots containing the mark divided by 
the average intensity of all spots not containing the mark. For 
peptide pull-down assays, 1 |jLg of recombinant proteins was 
incubated with 1 |jLg of biotinylated H3 or H4 histone peptides 
(Millipore, catalog nos. 12-403 for H3, 12-402 for Ac-H3, 12-372 
for H4, and 12-379 for Ac-H4) in 1 X HC buffer containing 200 mM 
KCl overnight at 4°C and washed extensively. Bound fractions 
were resolved by SDS-PAGE and visualized by silver staining. 

Yeast two-hybrid assay 

Full-length Bdf2, Bdf2-AN, Bdf2-AC, and Bdf2-BD were cloned 
into the Xmal/BamHI site of pGBT9 (Clontech) to generate 
fusion with the GAL4 DNA-binding domain. Epel was cloned 
into the Pstl/Bglll site of pGAD424 (Clontech) to generate fusion 
with the GAL4 activation domain. Both plasmids were trans- 
formed into the budding yeast strain pJ69-4A, and transformants 
were selected on medium lacking tryptophan and leucine to 
maintain both plasmids. The interaction of the two proteins was 
indicated by the activation of a HISS reporter, allowing growth 
on medium lacking histidine. 

RNA extraction and RT-PCR 

Total cellular RNA was isolated from log-phase cells using 
MasterPure yeast RNA purification kit (Epicentre) according to 
the manufacturer's protocol. Quantification with real-time RT- 
PCR was performed with Power SYBR Green RNA-to-CT one- 
step kit (Applied Biosystems). RNA serial dilutions were used as 
a template to generate a standard curve of amplification for each 



pair of primers, and the relative concentration of the target 
sequence was calculated accordingly. An actl fragment served as 
a reference to normalize the concentration of samples. The con- 
centration of each target gene in wild type was arbitrarily set to 1 
and served as a reference for other samples. 

Sir2 HDAC assay 

Full-length Sir2 cDNA was cloned into the pRSET bacterial 
expression vector (Invitrogen). Recombinant His-tagged Sir2 
protein was purified with Talon Metal Affinity Resin (Clontech) 
according to the manufacturer's protocol. Sir2 HDAC assay was 
performed using the SIRTainty Class III HDAC assay kit (Milli- 
pore) according to the manufacturer's protocol. Briefly, 5 [ig of 
recombinant Bdf2 or Bdf2-2YA proteins was incubated with 
P-NAD, nicotinamidase, tetra-acetylated H4 peptides, and dif- 
ferent amounts of recombinant Sir2 for 1 h at 37°C followed by 
further incubation with developer reagent for 1 h at room 
temperature. Fluorescence intensity was measured with a BioTek 
Synergy 4 Hybrid Microplate Reader with filter set to excitation 
at 420 nm and emission at 470 nm. 
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